Nonlinear Dynamics of Incoherent Superstrong Radiation in a Plasma 
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We present a new concept of nonlinear dynamics of incoherent superstrong radiation in plasmas. 
Recently we have disclosed a novel mechanism of the establishment of equilibrium between a photon 
and a dense photon bunch through the exchange of longitudinal photons (Tsintsadze 2004 Phys. 
Plasmas 11 855). Based on this mechanism of the "Compton" scattering type, we have generalized 
Wigner-Moyal equation for the dense photon gas, including the collision integral for the occupation 
number of photons. In the geometric optics approximation the Wigner-Moyal type of equation 
reduces to the one particle Vlasov-Boltzmann equation for the photon gas. From this equation, 
which gives a microscopic description of the photon gas, we derive a set of fluid equations, and 
consider numerically the formation of 3D shock waves. 



I. INTRODUCTION 



A general new approach to the study of the interaction of spectrally broad and relativistically intense electromagnetic 
(EM) wave packets (photons) with a matter was developed in our previous papers [1], [2]. This was inspired by the 
fact that the bandwidth of an initially coherent spectrum may eventually broaden, as a result of the several kinds of 
instabilities and the nonlinear interactions, and naturally for ultrashort pulses the initial bandwidth is increasingly 
broad. Moreover, in astronomical plasmas there are a variety of sources of radiation, and in this case we speak about 
the average density of radiation from all the sources and their spectral distribution. It should be emphasized that the 
radiation can be in two distinct states. Namely, one is when the total number of photons is not conserved. A good 
^v^j ' example is a black-body radiation. Another situation is when the total number of photons is conserved. Both these 
states have been studied in several aspects, mostly for the weak radiation. In spite of the abundance of publications in 
["~-. this field, a class of problems remains highly complicated and speculative, particularly in the case when the intensity 
t-H ' of radiation increases. Thus the study of property of such radiation (high-power, short pulse lasers with intensities 
up to l0 21 Wcm~ 2 , in which the photon density is of the order of 10 30 cm -3 , non-thermal equilibrium cosmic field 
radiation, etc.) is of vital importance. As was shown in Ref. [3] the photons acquire the rest mass and become one 
of the Bosons in plasmas and posses all characteristics of nonzero rest mass, i.e., we may say that the photon is the 
■ elementary particle of the optical field. 

In the past a several reviews have been published on theory of Bose-Einstein condensation (BEC) in a quantum Bose 
liquid [4] and in trapped gases [5]. Kompaneets[6] has demonstrated that the equilibrium between photons and electrons 
can be established via the Compton effect. Zel'dovich and Levich[7], using the kinetic equation of Kompaneets, have 
shown that in the absence of absorption the photons undergo BEC. Recently it was disclosed that exists an another 
mechanism of the creation of equilibrium state and BEC in a nonideal dense photon gas [8]. More importantly a 
new effect was predicted in the same paper, namely that the inhomogeneous dense photon gas can be found in the 
intermediate state. In the present paper, we consider the BEC and evaporation of the transverse photons (photons) 
from the Bose condensate. In our study we assume that the intensity of radiation is sufficiently large, so that the 
photon-photon interaction is more likely than the photon-electron interaction. In the case when the variation of the 
plasma density can be neglected in comparison with the variation of the photon density, the elementary excitations 
represent the longitudinal photons (photonikos), for which we will derive the well known Bogoliubov's energy spectrum. 



II. BASIC EQUATIONS AND STABILITY OF PHOTON FLOW 

We have shown in Ref. [3] that the behavior of photons in a plasma is radically different from the one in a vacuum. 
Namely, plasma particles perform oscillatory motion in the field of EM waves affecting the radiation field. The 
oscillation of electrons in an isotropic homogeneous plasma leads to the index of refraction, which depends on the 
frequency of the radiation, and is not close to unity for a dense plasma, i.e. 



where u> p = w A ' K< ^°' for an electron-ion plasma (neglecting the ion contribution) and uj p = «/ 2 • 4:7r ^" " for an 

electron-positron plasma (— e,mo e , riQ e and 7 are electron charge, the rest mass, density and the relativistic gamma 
factor of the electrons respectively). 
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Rewriting Eq.(l) in terms of an energy e — hu> and momentum p — hk (where h is the Planck constant divided by 
2ir) and introducing m 7 = Hlo p /c 2 , we obtain the expression of energy for a single photon 



e 7 = c(p 2 + rn 7 c 2 ^ = m 7 c 2 - (2) 



V2 „/ U 2 \-l/2 

c' 

which is expressed through the standard formula for the energy transport velocity 



M 7 = C| 



For the momentum of photon we can write 



Pi 



— hk — m 7 7 7 M 7 = m 7 ^l u 7 . (4) 



The form of Eq.(2) coincides with the expression for the total relativistic energy of massive particles, so that the rest 
mass m 7 is associated with the photon in a plasma. We note here that in contrast to material particles, the rest mass 
of photons depends on the plasma density. In view of this analogy between a photon in a plasma and a free material 
particle, we can treat the photon gas being with the plasma as a neutral subsystem of particles that have non-zero 
rest mass. 

In the following we will show that under certain conditions the photon-photon interactions dominate the photon- 
plasma particle interactions. In this case, we may consider the medium to consist of two weakly interacting subsystems: 
the photon gas and the plasma, which slowly exchange energy between each other. In photon-photon interaction the 
phases of the waves are, in general, random functions of time. Hence, we can average over them. In such a situation 
the perturbation state of the photon gas can be described in terms of the occupation number N(k, u), r, t) of photons, 
which is the slowly varying function in space and time. Recently, a new version of kinetic equation for the occupation 
number of photons, for modes propagating with the wavevector k and the frequency u), at the position r and time t, 
was derived in Rcfs.[l],[2],[9]. This equation reads 

| AT(fc, W) f, t) + ^(k- V)N(k, w, r, t) ^sin- (v ? ■ - | A) ■ p N ^J^ = , (5) 

where p = n e and n 0e are the non-equilibrium and equilibrium densities of the electrons, respectively, and 7 

can be expressed as 



r / n f dk f duj N(k,Lo,f,t) 

7=^=^1 + ^/^/^^, (6) 

where 3 — — 2hu>p 2 . 

Equation (5) represents the generalization of the Wigncr-Moyal equation (WME) for the classical EM field. It should 
be emphasized that from this equation follows conservation of the total number of photons, but not the momentum 
and energy of photons, i.e., 

N = 2 J dr J J ^-N(k,to,r,t) = const. , (7) 

where coefficient 2 denotes two possible polarization of the photons. Hence, the chemical potential, /i 7 , of the photon 
gas is not zero. 

If the local frequency to and wavevector k are related by the dispersion equation, u> = uj{k), then the occupation 
number is represented as N(k,u,r,t) = 2TrN(k,r,i)8(u> — to(k)) and Eq.(5) reduces to 

§-N(k, f, t) + C 2 W^M _ ^ 1 (v? . v s ) • = . (8) 

In the geometric optics approximation (sinx w x), Eq.(5) reduces to the one-particle Liouville-Vlasov equation with 
an additional term 
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It is important to emphasize that in Eqs.(5), (8), (9) there are two forces of distinct nature, which can change the 
occupation number of photons. One force appears due to the redistribution of electrons in space, Vn e , and time, 
^f-. The other force arises by the variation of the shape of wavepacket. In other words, this force originates from 
the alteration of the average kinetic energy of the electron oscillating in a rapidly varying field of EM waves, and is 
proportional to and J^i. Clearly, in the case when the density of plasma does not (or almost does not) change, 
remains only the force due to the variation of the shape of wavepacket. So that in the relativistic case the variation 
of the intensity of EM waves leads to the effects as self-modulation, self-focusing, formation of shock waves, etc. 

We now consider the propagation of small perturbations in a homogeneous photon flux-plasma medium assuming 
the temperature of electrons to be nonrelativistic, i.e., T e <C m 0e c 2 7. We linearize equations of motion and continuity 
of electrons with respect to perturbations, and seek plane wave solutions proportional to expi(qr — fit). Here we 
consider the range of low frequencies for which the inequalities f2 < qc < uo pe are fulfilled. In this case the photon 
flow can no longer excite the Langmiur plasma waves, and the contribution of perturbation of the electron density is 
rather small in comparison to perturbation of the photon density as 

Therefore, in our consideration the density of the plasma particles remain constant. In this case Eqs. (5) and (8) 
contain the variation of the relativistic 7 factor alone. 

We next demonstrate a new phenomena in the photon gas, which originates from the variation of shape of the 
wavepacket. Obviously, any waves, which can arise in the photon gas-plasma medium in this case, are the proper 
waves with the weak decrement of the photon gas. To this end, we linearize Eq.(8) with respect to the perturbation, 
which is represented as N(k,u(k),r,t) = No(k,u)(k)) + 5N(k,u(k),r,t) . The result is 

{— - n ) SN - 2_, (2i + l)!2«+i W (fc) ' 



or after summation we obtain 



_ n)§N= _ u2 S,<K(k + g/2) _ N-jk-g/2), 
V w / P 7o 1 u{k + q/2) u(k-q/2) > 

where £7 = = ^ I ' ~~~^frp~~^- After integration over all wavevectors k, from Eq.(12) we get the 

dispersion relation due to the relativistic selfmodulation 

, + ± f f^lWfim _ £M?1 W (| - n) _ , (13) 

27o J 1 u(k + q/2) u(k-q/2) > w Vw(fc) / 

which has, in general, complex roots. In Eq.(13) N (k) is the occupation number in the equilibrium state and is 
represented as 



N (k) = no 7 (2^ ) ~ 3/2 ^(-^^) • (14) 



This is a spectral Gaussian distribution, with the average wavevector k and the spectral width at- 

We can rewrite Eq.(13) in another form, taking into account a pole in the integral — qu = , where u — < c. 
Using the well known relation 

lim — - — = p— — nr5(x) , (15) 

where p denotes the prescription that at the singularity x — the principal value is to be taken, Eq.(13) is rewritten 
in the form 

i + ffi p I + m - N °~$ = m }/(** - ") 



27o J u{k)\ uj{k + q/2) ' u(k-q/2) 
_ J r d 3 k ( N+(k + q/2) _ Ny(k-g/2) } 
27^ J u(k)\ uj(k + q/2) w(k-q/2) > 



+ ^ J <PkrN+(k + g/2) _ N -(k-g/2) ]m _ = Q _ 
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In the first and the third integrals we now replace the wavevector k + q/2 by k, and in the second and the forth 
integrals k — q/2 — * k. Assuming that k = k + x, \ q \, \ X 1^1 1 7 after integration we obtain the dispersion relation 

i + fXi ^^ ( diuxA ) 

+ (n - qu g ? - Q 2 V S 2 - a 2 q* + 2 7o 3 c 2 u;(fco) V d Xz ) x , = ^(fi-*8 9 ) ' 1 ' 

where u g = qu g = qu g cosO, N ( Xz ) = J d Xx J d Xy N (x), V E = c^^y^i, a = and the 

velocity V s — cv^^^y can be treated as the sound velocity of longitudinal photons, similar to phonons in quantum 
liquid at almost zero temperature [4]. 

We now neglect the small imaginary term in Eq.(17), and examine the following equation 

{ n~qu g ) 2 =q 2 {v 2 + a 2 q 2 -V 2 }. (18) 

This equation has a complex solutions, when the inequality > V 2 + a 2 q 2 , for q < 2u} ^ K 2 > holds. The 

unstable solution can be written as £1 = £1/ + iCl//, where 

il/ = qu g cos0 (19) 

and the growth rate for the unstable modes is 

£1// = qy/v% - V 2 - a 2 q 2 . (20) 

This instability at V s = a = was disclosed in Ref.[2]. 

If £1/ >| ft// I, then the solution (19) clearly describes the emission of photonikos by the bunch of photons inside a 
resonance cone (cosO = Sl//qu g ), similar to the well known Cherenkov emission of EM waves by a charged particles. 
In addition, here is an analogy with the Landau criterion about the creation of elementary excitations. 

We now consider the case, when V 2 > Vj or u 2 (-fQ — 1) < &^l<y 2 c 2 . If we write Eq.(18) in the moving frame 
(ft — qu — > ft), and multiply the both sides by h, then we obtain the well known Bogoliubov energy spectrum 



^=\h 2+ (2^Y> ^ 



where e — hCt, p — hq, and m e ff — froj/c 2 . This equation has the same form as the energy of the elementary 
excitations in the quantum liquid at zero temp erature. 

For small momentum (p < m e ffV — meff\JV 2 — V£.) from Eq.(21) follows e = pV w V s p. This expression reveals 
that the coefficient V is the velocity of "sound" in the photon gas, similar to phonons in the quantum medium. For the 
large momentum (p ^> m e ffV), the photoniko energy (21) tends to p 2 /2m e //, i.e., the kinetic energy of an individual 
photoniko. 

The decrement, Imfl = ft//, of photoniko is derived from Eq.(17), which is 



where Vq — c 2 J-s^-p F — yJ~^(^ J^k ) ) (3^) ^ e n0 * e ^ e decrement has a maximum near the frequency of 
the Cherenkov resonance O ~ qu cos Q. 

In the above we have derived the condition (10), which implies that the density of plasma almost does not change. 
In this case from Eq.(5) follows that photons with different frequencies and wavevectors scatter on the photon bunch, 
and the equilibrium state may be established. This is a new phenomena of the "Compton" scattering type, i.e., the 
statistical equilibrium between the photons and the photon bunch will establish itself as a result of the scattering 
processes. To better exhibit this, we derive the Pauli equation in the wavevector representation. Note that the 
Wigner-Moyal-Tsintsadze equation (WMTE) (5) is convenient for a study of non-linear processes in the photon gas 
in the case of so-called weak turbulence states. Such states of the photon gas-plasma means that the energy of the 
quasi-particle (photoniko) is small compared to the photon energies. In the case of weak-turbulence, we can define the 
distribution function N(k,u,r,t) in any approximation for the amplitude of perturbation by iterating the WMTE. 
The result is 

dN(k) ^ f d 3 q 



dt 



£ / ^[w±(k + q,k)N(k + fi-W±(k,k + q)N(k)] . (23) 
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Here we have introduced the scattering rate 

7T w£ I Sp(q) | 2 ^ 
W± = 77 — — * (H±g-«/) ■ (24) 

Equation (23) describes the three wave interaction. Namely, the photon passing through the photon bunch absorbs 
and emits photonikos, with frequencies f2 = ^(lu — loi) and wavevectors q = =p(fc — fc/). Integral in Eq.(23) is the 
elastic collision integral and describes the photon scattering process on the variation of shape of the photon bunch. 

It should be emphasized that the kinetic equation (23) is a new version of the Pauli equation in the wavevector 
representation [10]. Note that the equation type of (23) has been obtained by Pauli for a quantum system, whereas 
Eq.(23), derived for the dense photon gas, is pure classical. This equation indicates that the equilibrium of the photon 
gas is triggered by the perturbation Sp = 6(n/n , y), in particular by the perturbation in the wavepacket. 

III. BOSE-EINSTEIN CONDENSATION IN THE PHOTON GAS 

It is well known from plasma physics that the transverse photon can decay into the transverse photon and a plasmon 
(electron Langmuir and ion sound plasmons). For this process to take place it is necessary to have a fluctuation of the 
plasma density. In our consideration, as was shown in Sec. 2, the density of plasma remains constant. So that only 
one mechanism (relativistic effect), which can cause BEC, is the decay of the transverse photon into the transverse and 
the longitudinal photons. Namely, the photon (u>, k) generates photoniko and scatters on latter, with the frequency ui 
and the wavevector kl. In this case, the wave number and frequency matching conditions are satisfied, i.e., k = kl + q, 
and uj = wl + O. These processes continue as a cascade (uoi = u)ll + ft/, kl = kn + ql, etc.) till the wavevector of 
the photon becomes zero and the frequency u> = ujp/"/ 1 ^ 2 . This cascade leads to the BEC of the photons and the 
creation of the photoniko gas. When the density of latter is sufficiently small, the photonikos may be regarded as 
non-interacting with each other, however they can exchange energies and momentum with a bunch of the photons by 
scattering mechanism, and finally such process can lead to the equilibrium state of the photoniko gas, with the Bose 
distribution 

n L = , (25) 

where f2 is the Bogoliubov spectrum (21). We note here that since the total number of photonikos does not conserve, 
the chemical potential of this gas is zero. The result of BEC is that the ground state is filled by photons with the total 
rest energy E = N ^y^, where N is the total number of photons in the ground state. We can express explicitly the 
energy of the ground state, Eq, through the total number and the volume of the photon gas in two limits. First, we 
consider ultrarelativistic case, when 7 ss Q 1 / 2 = (hcJ p /m c 2 ) 2 / 3 (N /N 0e ) 2 ^ 3 (l/V 1 / 3 ), where uJ p = (4we 2 N 0e /m 0e ) 1 / 2 , 
-/V 0e is the total number of electrons. 

The chemical potential of the photon gas is 

Note that in quantum Bose liquid at T = 0, the chemical potential increases with increase of A*o, whereas in the 
present case p decreases as shows Eq.(26). 

We now find the pressure in the condensate as P c = —dE /dV, and write the equation of state for given N , the 
result is 

P C V A/3 = const . (27) 

Next in the nonrclativistic limit, we have 7 ~ 1 + (l/2)(hLJ p /m c 2 )(N /N 0e )(l/V 1 / 2 ). Then for the equation of state 
and the chemical potential we get the following expressions 

P C V 3/2 = const , (28) 



U = TlUJ p ( 1 - . 29 
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Note that the second term in Eq.(29) is always less than one. 

The kinetic equation (23) is the simplest nonequilibrium model in which we can expect BEC. To show this, we 
introduce the number density of photons in the condensate as n = limy^oo Then we can write N(k,t) = 
4ir 3 S(k)no(t), which means that the photons with a zero wavevector are in the condensate. 

We now discuss an another situation in which the photon with the wavevector k is collided with the photoniko with 
q. For such processes exists the probability of the condition k + q = to be satisfied. Thus the photons absorbing 
the photoniko will pass to the ground state. Therefore, the occupation number, which is under integral in Eq.(23), 
can be written as 

N(k + q) = 4n 3 5(k + q) ni (t) + N(k + q)^_j , (30) 

where n\{t) is the density of pairs, for which the equality k + q = holds. Substituting Eq.(30) into Eq.(23), for the 
Bose-Einstein condensate density n a (t) we obtain 

dn (t) n 1 (t)-n (t) +2 f Sk, f A ^ _ ^ (31) 



dt t q J^oi^rJ (2tt) 3 

where Tq" 1 = / -0^W{q) = § lu* J ^-^0^s(n - ^j7 1/2 ) , and in the last term we have replaced k + q^kf. 

We specifically note here that the first and the third terms in Eq.(31) are the sources of production of the Bose- 
Einstein condensate, whereas the second term describes the evaporation of photons from the ground state. It is clear 
from Eq.(31) that the stationary solution may exist, and it is 

no = ni H n 2 , (32) 

< T > 

where < r > *= £ / with r{k,)^ = J $frW{kt, k, - q). 

The problem of BEC and evaporation of the Bosc-Einstcin condensate we can investigate by Fokker-Planck equation 
(FPE), which we derive from the Pauli equation. Here we suppose that the wavevector q and the frequency £1 of 
photonikos are small in comparison with the k and uo. This allows us to use the following expansion in the integrand 

- d a a d 2 WN 

W(k + q,k)N(k + q)^W(k)N(k) + q-^(WN) |, =0 U=o +■■■ ■ (33) 

Substituting this expression into Eq.(23), we obtain the FPE for photons, which describes the slow change of the 
occupation number in the wavevector space 

where 

d 3 q 



{ . )7 _ )3 q t W(q,k) , (35) 

D l0 = J -^ m3 W{q,k) . (36) 

The quantity A t is the dynamic friction coefficient, whereas D l} is the diffusion tensor in the wavevector space. 
Introducing the following definition B t = A t + ^§jf-, and noting that the expression on RHS in Eq.(34) is divergent 
in wavevector space, Eq.(34) can be rewritten as 

dN d 

where ] l — —B^ — -rpfj^- is the photon flux density in wavevector space. The fact that the flux should be zero 
allowed us to express B t and D ZJ in terms of one another. The equilibrium distribution function we choose to be 
Gaussian 

_ fc 2 

N(k) = const -e^l . (38) 
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Substituting expression (38) into the equation j — 0, we obtain B l a\ — ^kjD tJ . Finally the transport equation of the 
photon gas takes the form 



dt 



In order to solve Eq.(39), we consider a simple model representing the diffusion tensor as D v — DoS l0 , with Do = const. 
Such situation is realized when u ■ cos0 <C qc 2 /ui. With this assumption Eq.(39) reduces to 

f-|<»> + ^, (40, 

where a = We first neglect the diffusion term in Eq.(40), and assuming fc(0, 0, fc), we get for / = kN 

The general solution of this equation is an arbitrary function of fco = ke at (where fco is the initial value of the 
wavevector), i.e., /(fco) = f(ke at ). This function is constant on the curves fco = ke at , while at t — > oo the wavevector 
goes to zero (fc — > 0). The meaning of this is that the occupation number would tend to peak toward the origin as 

N(k,t) = {- = const -e at . (42) 
fc 

Thus we conclude that the friction effect leads to BEC of the photon gas. Here the condensate formation time, t c , is 
defined by the relation at c <~ 1. It should be emphasized that the exponential decay of the wavevector k = k e~ at is 
due to the linearity of the FPE. The presence of nonlinearity could saturate the exponential growth of N(k,t). 
Second, we suppose that W(q, fc) is the even function with respect to q. In this case A = 0, and we have 

Assuming that initially all photons are in the ground state with k = 0, and the occupation number of the photons is 
N = 47r 3 n (5(fc), then the solution of Eq.(43) reads 

N(k, t) = n a e~^ /(2irD t) 1/2 . (44) 

In the course of time, the number of photons with fc = decreases as i~ 1//2 . The number of photons in the surrounding 
wavevector space rises correspondingly, and initially peaked at the origin is to be flatten out. Let us determine the 
mean square wavevector from origin at time t. From expression (44) we have 

< fc 2 >= D Q t . (45) 



Thus, V < fc 2 > increases as the square root of time. Eqs. (44) and (45) manifest that the evaporation of photons 
from the condensate takes place. 

We now derive a relation between the diffusion time, to, and the time of the condensation. From Eqs. (42) and (44) 
we obtain 

* 2a l ; * fc2 



Do 2Dn 



where ct/c = as follows from the Wigncr function N(k,r) = Noexp(^—^pz — ^rj- Note that r is the initial cross 
size of the photon bunch. Finally, we arrive at the desired relation 

^=fc 2 r 2 . (46) 

From here it is evident that for the condensation and evaporation of photons, it is necessary that the following 
inequality t^ ^> t c is satisfied. This is to be expected, as in realistic astrophysical objects with a strong radiation, as 
well as in laboratory experiments with a laser radiation, the following condition fc 2 r 2 3> 1 is valid. 
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IV. FORMATION OF SHOCK WAVES 



We now consider 3D problem of the nonlinear stage of evolution of the photon bunch, revealing the shock formation. 
From the kinetic equation (9) one can obtain a set of fluid equations by the usual way. These equations of motion 
and Q of the photon gas are: 

^ + (S 7 .V)S 7 = -vi, (47) 



+ divQu 7 = , (48) 

where u 7 is the photon mean velocity. 

We have carried out 3D numerical simulation that self-consistently solves the above fluid equations. We first discuss 
numerical analysis for pancake-shaped photon beams, i.e., xq = yo = 2, and zq = 0.5. The propagation direction is 
taken to be in the z-dircction, and Q = 8. Figures la and lb show that the width of the beam becomes shorter 
in time along the propagation direction. The photons are concentrated in the small region (on the z-axis), and at 
t=1.8 starts to form the shock wave. Whereas, Figures lc and Id show the insignificant expansion of the beam in the 
transverse direction. 

Rather different processes develop when the initially the photon beam has a form of the bullet, i.e., when xq = 
yo = 0.5, z = 2, and Q — 8. Namely, as can be seen from Figures 2a and 2b, the beam elongates in time along the 
propagation direction. Whereas, the compression of the beam in the transverse direction is observed in Figures 2c 
and 2d, and the shock wave takes place at t=1.6, in a time short compared to the previous case. 



V. SUMMARY AND DISCUSSIONS 



We have investigated the interaction of spectrally broad and relativistically intense EM radiation with a plasma 
in the case when the photon-photon interaction dominates the photon-particle interactions. We have established the 
condition under which the variation of the plasma density can be neglected in comparison with the variation of the 
photon density. In such case the elementary excitations represent the photonikos, for which we have derived the 
well known Bogoliubov energy spectrum. We have studied the BEC and evaporation of the photons from the Bosc 
condensate. To this end, from the WMTE we have derived the Pauli kinetic equation for the photon gas. We have 
also derived the FPE, and presented a simple model, which exhibits the possibility of the creation of Bose-Einstein 
condensate and evaporation of photons from the condensate. The theory developed in this paper is of essential interest 
to study an extremely complex phenomena of astrophysics in laboratories using ultra-intense lasers, ft should also be 
the case for a black hole. In this connection we speculate that the recently observed radiation from the black hole may 
be attributed to the evaporated photons from the Bosc condensate. That is initially all photons fall into Bose-Einstein 
condensate, and then after a certain time, as discussed in this paper, some photons undergo evaporation from the 
condensate. 
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